Richardson-Gaudin integrability in the contraction limit of the quasi-spin 
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Background: The reduced, level-independent, Bardeen-Cooper-Schrieffer Hamiltonian is exactly diagonalizable 
by means of a Bethe Ansatz wavefunction, provided the free variables in the Ansatz are the solutions of the set 
of Richardson-Gaudin equations. On the one side, the Bethe Ansatz is a simple product state of generalised 
pair operators. On the other hand, the Richardson-Gaudin equations are strongly coupled in a non-linear way, 
making them prone to singularities. Unfortunately, it is non-trivial to give a clear physical interpretation to the 
Richardson-Gaudin variables because no physical operator is directly related to the individual variables. 

Purpose: The purpose of this paper is to shed more light on the critical behavior of the Richardson-Gaudin 
equations, and how this is related to the product wave structure of the Bethe Ansatz. 

Method: A pseudo-deformation of the quasi-spin algebra is introduced, leading towards a Heisenberg-Weyl 
algebra in the contraction limit of the deformation parameter. This enables an adiabatic connection of the exact 
Bethe Ansatz eigenstates with pure bosonic multiphonon states. The physical interpretation of this approach is 
an adiabatic suppression of the Pauli exclusion principle. 

Results: The method is applied to a so-called "picket-fence" model for the BCS Hamiltonian, displaying a typical 
critical behavior in the Richardson-Gaudin variables. It was observed that the associated bosonic multiphonon 
states change collective nature at the critical interaction strengths of the Richardson-Gaudin equations. 

Conclusions: The Pauli exclusion principle is the main responsible for the critical behavior of the Richardson- 
Gaudin equations, which can be suppressed by means of a pseudo deformation of the quasispin algebra. 

PACS numbers: 02. 30. Ik, 21.10.Re, 21.60.Ce, 74.20.Fg 



I. INTRODUCTION 

The description of superfluid properties in many- 
article systems such as metals [J Q or atomic nuclei 
0] involves the process of pairing. The pairing inter- 
action will cause two particles in time-reversed single- 
article states to gain energy by forming a Cooper pair 
In the thermodynamic limit, the ground state of 
the system will constitute a condensate of Cooper pairs, 
which is well approximated by the Bardeen-Cooper- 
Schrieffer (BCS) coherent state description of supercon- 
ductivity [l|, Q . While the approximation is tailor made 
for the thermodynamic limit, it is less appropriate in the 
mesoscopic finite-size regime @,@]- The major source for 
inaccuracies can be found in the particle-number fluctu- 
ations in the BCS state, which become relatively strong 
with respect to the total number of active particles. 
Thus, the ground-state structure of finite-size systems is 
more intricate than a straightforward BCS condensation 
of Cooper pairs, so there is a clear call for canonical ap- 
proaches. Direct numerical diagonalisation algorithms, 
such as e.g. Lanczos, can easily deal with small systems, 
but they play a loosing game for large finite-size systems, 
as is the case in e.g. superconducting metallic nanograms 
or exotic atomic nuclei Q , due to the exponential size 
explosion of the Hilbert space. 
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Remarkably, the reduced BCS Hamiltonian, with level- 
independent scattering terms, was found to be integrable 
fiol . A common characteristic of integrable systems is 
that they support as many (non-trivial) conserved oper- 
ators, commuting with the Hamiltonian, as there are de- 
grees of freedom in the system Earlier, Richard- 
son had shown that the reduced BCS Hamiltonian can be 
diagonalised exactly by means of a Bethe Ansatz wave- 
function, provided the free variables in the Ansatz (also 
referred to as rapidities [l4[ or pairons [Hj]) form a so- 
lution to the set of non- linear coupled equations, the so- 
called Richardson-Gaudin (RG) equations Once 
the RG variables have been determined, all other physical 
observables can be evaluated from these. 

Within the framework of the Algebraic Bethe Ansatz, 
the eigenstates of an integrable model are given by prod- 
ucts of generalised ladder operators [l?} • Whereas these 
generalised operators can become quite involved in terms 
of the physical operators acting on the Hilbert space 
fl8| . the Bethe Ansatz eigenstates of the reduced BCS 
Hamiltonian are product states of generalised pair opera- 
tors. The factorization of Richardson's Bethe Ansatz into 
physical operators instigated several investigations re- 
garding its relation to BCS Cooper pairs projected- 
BCS Cooper pairs @, H|, and bosonic excitations [2fjj |. 
The clear correspondence between the latter structures 
and Richardson's exact solution is often obscured by the 
interpretation of the RG variables. For one reason, it 
is well-established that the RG equations display a crit- 
ical behavior at certain values of the pairing interaction 
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strength, although the energy eigenvalues of the Hamil- 
tonian behave analytically |2l| - t23| . Also, no physical ob- 
servable has been identified yet that can probe the indi- 
vidual RG variables as physical entities, nevertheless, it 
has been noticed that a qualitative interpretation of the 
RG variables in terms of the collective behavior of the 
associated pair creation operators is possible @, H(| ■ 

The purpose of the present paper is to shed more 
light on the critical and collective behavior of the RG 
variables by making an adiabatic connection with pure- 
bosonic product state. This will be done by means of 
a pseudo-deformation of the su(2) quasispin algebra, re- 
lated to a suppression of the Pauli exclusion principle. 
The intcgrability of the reduced BCS Hamiltonian and 
the Richardson- Gaudin formalism will be briefly pre- 
sented in the next section (Sec. |TTJ) . The following sec- 
tion is divided into four parts, containing a recapitula- 
tion of the required ingredients of the pp-T&mm Dancoff 
Approximation (Sec. MI A]) , a discussion of the algebraic 
properties of the pseudo-deformed algebra (Sec. IIIIBD . 
its embedding into Richardson-Gaudin intcgrability (Sec. 
IIII Cj) . and how we can understand more about the crit- 
ical behavior of the Richardson-Gaudin equations (Sec. 
IIII Dl) . Conclusions are presented in the final section (Sec. 

ITB. 



II. RICHARDSON-GAUDIN 

The reduced BCS Hamiltonian is given by 

rn m 

H = J2^i + gJ2 Sj^, (l) 

i— 1 ij — 1 

where the Roman indices {i, j = 1 . . . m} are shorthand 
for a set of single-particle quantum numbers, with single- 
particle energy e, and degeneracy Hi. Within the spher- 
ical shell model, the degeneracy is Hi = 2ji + 1 where 
ji is the total angular momentum of the level i. When 
no spherical symmetry is implied, Hi = 2 refers to the 
two-fold degeneracy of the intrinsic spin. The Hamilto- 
nian ([TJ) consists of single-particle operators, counting the 
number of fermions within a level i, 

hi = 2J (oL^mi + a\- ni a m ), (2) 

mi>0 

and a fermion-pair scattering component, represented by 
the pair creation/annihilation operators 

Sj = a^.a^., Si — (Sj)^ = cimMjn^ (3) 

m^>0 m^>0 

with a^. (a mi ) the standard fermion creation (annihi- 
lation) operators. The bar notation denotes the 
time-reversed partner of the corresponding operator, e.g. 
= (^)"' i_mi a]. _ TO . in the spherical shell model @. 



The set of operators (JSJ)-© span an sit (2) quasispin al- 
gebra 

[S^ , Sj^ = SijS^j [S>i , Sj] — SijSi, [^i: $j] = %5ijSi , 

(4) 

where Sf = ^fii — k^i- The irreducible representations 
(irreps) are given by 

\di,fj,i) = \lCli- \vi,\m - iOj), (5) 

with Vi the number of unpaired fermions or seniority [3] , 
such that the total number of particles within a level i 
is given by rii = 2Ni + Vi with iVj the number of pairs. 
Richardson's result states that the reduced BCS Hamil- 
tonian ([TJ) can be diagonalised by means of a product 
state of generalised pairs acting on the 0™ x su(2)i low- 
est weight state \9) = YiiLi \di,—di) 

N / m Af \ 

l« = II Edr (8) 

Q = l \i = l 1 a J 

provided the set of RG variables {x} form the solution of 
the RG equations 

m , N 1 

1 + 2,^-^-2.9 Y, ( 7 ) 

i=i £i Xa 0=1,?* Xf3 a 

for a = 1, . . . , N with N the total number of pairs. The 
curly bracket notation {x} = {xi, x%, . . . , Xjsr} will be 
used to denote a set of RG variables, identifying a state. 
Once the RG variables {x} have been determined, the 
total energy of the eigenstate is given by 

N m 

E = ^T / x a +^2e l v l . (8) 

a=l i=l 

There are multiple strategies to derive this result, ei- 
ther via Richardson's original approach pHHH, Gaudin's 
integrability constraints jToj . a linearization of the R- 
matrix in the Algebraic Bethe Ansatz [23, or via a com- 
mutator scheme using a Gaudin algebra j(| [26| . A brief 
discussion of the latter approach can be found in Ap- 
pendix [A] in support of the results of Sec. IIII CI 

III. HARD-CORE, GENUINE, AND BOSONS IN 
BETWEEN 

A. The Tamm-Dancoff Approximation 

As already mentioned, Richardson's Bethe Ansatz 
state ([5]) is a simple product state of generalised pair- 
creation operators, which are highly correlated in a non- 
linear way by means of the RG equations ([7]) . A sensible 
approximation would be to relax the intricate correla- 
tions described by the RG equations by a simpler and 
more tractable prescription, retaining the product struc- 
ture of the wavefunction. This is exactly what is achieved 
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by certain factorization approximations, such as the pro- 
jected BCS [27|, the Random-Phase Approximation [28[ 
or the pp- Tamm-Dancoff Approximation (pp-TDA) (cfr. 
Chap. 8.2.3 of H3] ) . The basic assumption of the pp-TDA 
is that the eigenstates of a Hamiltonian can be consid- 
ered as harmonic excitations (multiphonon states) of el- 
ementary eigen modes (phonons) of the system. For the 
pairing Hamiltonian, viable candidates for the elemen- 
tary excitation modes are given by the 1-pair excitations 



Hj2Y i Sl\9} = tuJ2 Y i8l\e)- (9) 

1=1 i=l 

In the case of the reduced BCS Hamiltonian ([1]), the 
eigenvalue equation ([9]) has a geometric interpretation 
because the m eigenmodes hu>k (k — l,...,m) can be 
determined by the zeros of the secular equation of the 
pp-TDA 

i + ^^ = °. do) 

i— 1 

The coefficients of the fcth eigenmode are given by 
Y* = 2e --hu ' so ^ ne -^"P a i r wavefunctions can be read- 
ily constructed 

JV / m A t \ 

where {huj} is shorthand for a set of pair variables 
{fadMu, . . . ,hu}f.(w)} with k(a) an integer- valued func- 
tion selecting the fcth eigenmode for the cvth pair. In 
the coming, the set {hut}, defining the state (jlip will 
be represented by the TDA distribution (u\, 1/2, • ■ • , Vm) 
where t>k denotes the number of generalised pairs with 
eigenmode haj^ in the product state (fTTj). Based on the 
geometric interpretation of the secular equation, it is seen 
that the eigenmodes fwjk are bound between two subse- 
quent single-particle levels 2ek-i < Swfc < Sk, except 
for hw\. which is bound above by 2ei but not bound 
below (see e.g. Fig. 6.6 in 0). This eigenmode is com- 
monly referred to as the collective eigenmode, and a mul- 
tiphonon state (fTTj) with distribution (N, 0, . . . , 0) will 
also be called collective. 

Given the exact solution of the problem, it is straight- 
forward to check the fitness of the approximation, either 
by confronting the exact energy spectrum with the ex- 
pectation value of the Hamiltonian calculated with the 
multiphonon reference state (fTTj) . or by calculating the 
overlaps of this state with the exact Bethe Ansatz ([5]). 
The latter has been done by Sambataro [2(| , where a so- 
called "picket-fence" model [29| of 12 uniformly spaced 
two-fold degenerate single-particle levels has been con- 
sidered in which 12 paired fermions (N = 6) are ac- 
tive. The real and imaginary part of the RG variables 
are plotted in Fig. [TJ showing a typical behavior. All 
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FIG. 1. The real (a) and imaginary (b) part of the RG vari- 
ables of the "picket-fence" model, employed in [S^l, as a func- 
tion of the interaction strength \g\. The singular interaction 
strengths are highlighted by means of dashed lines. Dotted 
and dashed lines identify the interaction strengths where the 
TDA distribution in the contraction limit, denoted by the 
vectors (yi, V2, Vm), changes nature (see Sec. IIII D| . 



RG variables are real valued for sufficiently weak inter- 
action strength (\g\ < 0.393 in Fig. [TJ. As the interac- 
tion strength increases, the variables recombine two by 
two to form complex conjugate pairs at the critical in- 
teraction strengths g Cl until all variables have become 
complex conjugate pairs (\g\ > 0.891). The overlaps 
(tj}({x})\tp({hLij})} were calculated in [20| for the ground 
state and a few excited states as a function of the interac- 
tion strength g. The conclusion for the ground-state over- 
laps was that a distinction can be made between 3 differ- 
ent regimes: a weak-coupling (\g\ < 0.25), intermediate- 
coupling (0.25 < \g\ < 1.25) and a strong-coupling regime 
(1.25 < \g\). It was observed that the overlap with the 
ground state was largest with the (1, 1, 1, 1, 1, 1, 0, . . . , 0) 
multiphonon state for the weak-coupling regime whereas 
the overlap was largest with the collective (6,0, ...,0) 
multiphonon state in the strong-coupling regime. For 
the intermediate regime, the situation was more ambigu- 
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ous as no single multiphonon state was found to dom- 
inate over the others. To illustrate the situation, the 
overlaps of the ground state with a selected set of mul- 
tiphonon distributions are given in Fig. [2] The distri- 
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FIG. 2. Overlaps of the normalized ground state |^({a;})} 
(HU of the reduced BCS Hamiltonian {JXJ) with the normalized 
product states ^fffio;l')) (|ll[l for the model described in [20| . 
The notation (ii,i2, . . . ,i m ) denotes how many generalised 
pairs ik are chosen with frio^. 

butions have been chosen as such that the distribution 
with maximum overlap with the exact ground state is 
present in the figure for every interaction strength g. It 
is readily seen from the figure that several multiphonon 
states compete for the largest overlap in the intermediate- 
coupling regime, pointing out that the structure of the ex- 
act ground state is quite intricate, opposed to the weak- 
and strong-coupling regime where a clear picture of the 
exact ground state emerges, thanks to the dominant over- 
lap of a well-defined multiphonon state. It can be seen 
from Fig. [5] that the maximum overlap changes from 
(1,1,1,1,1,1,...) to (1,1,1,1,2,0,...) at \g\ = 0.385 in 
the weak-coupling regime and from (5, 1, 0, 0, 0, 0, ... ) to 
(6,0,0,0,0,0,...) at \g\ = 0.872 in the strong-coupling 
regime. These values are remarkably close to respectively 
the lowest and highest critical points of the RG equations 
(see Fig.[T|), so it would be interesting to know if the crit- 
ical point in the intermediate regime at \g\ = 0.603 also 
corresponds to a change in multiphonon nature of the pp- 
TDA. Unfortunately, the picture is quite blurred in the 
intermediate regime, so a different approach is needed to 
enable a clear connection. This will be discussed in the 
following subsections. 



B. A pseudo deformation of the quasi spin 

Fermionic pairs are similar to bosons in the sense that 
they mutually commute. However, opposed to bosons, 
fermionic pairs feel the presence of other pairs through 
the Pauli exclusion principle between the constituent 
fermions. For these reasons, fermion pairs are sometimes 



referred to as hard core bosons. It is possible to transform 
the hard core bosons into genuine bosons by introducing 
a pseudo deformation of the quasispin algebra. By re- 
defining the pair creation and annihilation operator such 
that the influence of the number operator is tunable by 
means of a deformation parameter £ € [0, 1] 



$(0.£i(0]=2fci(^ 



(12) 



we arrive at a pseudo deformation of the quasispin alge- 
bra, as introduced by Arecchi et. 



f the que 

mini 



[s}(0,sm] = M£2#(0 + (£ - 1)10*), (i3) 



with 



IQ. 



a ^-independent operator and 



hi retaining its interpretation as number operator. It is 
worth emphasizing that no underlying structure is as- 
sumed for the creation- and annihilation operators S\ (£) 
and 5*i(£), in contrary to the regular quasispin algebra 
(see Eq. [3|) , so these operators are solely defined by means 
of the algebra (|12p . The original quasispin algebra (fj| is 
recaptured in the £ = 1 limit of the algebra, and a non- 
normalized bosonic Heisenbcrg-Weyl algebra hw(l) is ob- 
tained in the £ = contraction limit [301 ] . It is preferable 
to employ the term pseudo deformation as the deformed 
algebra (fl~3|) effectively spans an su(2)^ algebra, defined 
by the generators 



Si(0 



m) = ^sj(o, mo 



(14) 
(15) 



Thanks to the explicit form of the deformed Cartan op- 
erator ^4°(£), it is possible to assign a loose physical in- 
terpretation to the deformation parameter £. The action 
of the Cartan operator on the lowest weight pair vacuum 
state \di(£),—di(£)) yields the following relation for the 
su(2)^ irrcp label 



dm 



4£ J L l 



(16) 



As a result, the deformation parameter £ opens up a given 
shell by increasing the effective degeneracy admit- 
ting more pairs than strictly allowed in the original irrep 



di(l) 



4 J 



In the contraction limit of £ — > 0, the 



effective degeneracy reaches infinity, which is consistent 
with a bosonic interpretation. It should be emphasized 
that the assignment of c?i(£) as an su(2)^ irrep label is 
only physically meaningful as long as the deformation 
parameter is £ = 1+9 \./ n . with k G N; otherwise, the 
irreps are no longer unitary [32j |. Nevertheless, this re- 
quirement can be relaxed in the framework of Richardson 
integrability, as everything is well-defined from the com- 
mutation relations (|12[) and the action of the operators on 
the pair vacuum for the deformation parameter over the 
regime £ e [0, 1]. Fig. [3] illustrates how the su(2)^ irrep 
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FIG. 3. An illustration of the sm(2)j quasispin irreps with 
label for Qi — 6 and Vi = 1. The unitary irreps, at £ = 

i+2fc/n- with fc € N, are denoted by open circles, connected 
by dashed lines. The lines connecting the lowest and highest 
weights with = T^(0 are drawn in full. For f = 1, 

the = 1 irrep is retained. The operators A\ and At, 

depicted in the su(2)^ root diagram, are defined in Eq. (|14| 



provided the set of deformed RG variables form 
the solution of the set of deformed RG equations 



2ff£ 



N 

E 



1 



Z/3 



0. (21) 



Remark that the value £<i(£) = j^i - £|^i i s well-defined 
and finite over the whole regime £ £ [0,1]. There are two 
ways to obtain this result. The first method is rather 
straightforward, as discussed in Appendix \K\ One sim- 
ply commutes the Hamiltonian (|19|) through the product 
state ([2U1) using the commutation relations (fl~2"|) (or cquiv- 
alcntly (|13p ). and requires the non-diagonal byproducts 
to vanish. The second approach employs the effective 
su(2) ? algebra dD for £ G]0,1]. The Hamiltonian (O 
can be rewritten as 

ra tci 

£(0 = 5>(2i?(0+£n)+s£ £ 4(0^(0, ( 22 ) 



label dj(£) grows with decreasing £, until the contraction 
limit, where the algebraic structure changes from su(2)^ 
to hw(l). 

The algebra (|14|) becomes ill-defined in the contraction 
limit £ = 0. For this limit, it is better to resort to the 
original algebra. By defining the normalized operators 



Si(0), 



S?(0), (17) 



the bosonic nature of the creation/annihilation operators 
in the contraction limit becomes apparent 



<%st, 



[s?,s. 



sti 



(18) 



The operators s? can be considered equivalent to the 
number operators 



si Si 



a-. In the fol- 



lowing subsection, the integrability of the reduced BCS 
Hamiltonian ([1} will be investigated using the present 
pseudo deformation of the quasispin algebra. 



C. The pseudo deformed RG equations 

The key observation is that the reduced BCS Hamilto- 
nian ([TJ remains intcgrable when the deformed commu- 
tation relations arc used. The Hamiltonian 



i—l i,j=l 

is diagonalizablc by means of the product state 

N / m ^f. 



(19) 



MM®})) = n e 



sm 



(20) 



which can be interpreted as a reduced BCS Hamiltonian, 
entirely expressed in terms of su(2)^ generators, with a 
scaled coupling constant gt; and deformed su(2)^ irrep la- 
bels HU). Filling out these values in the regular RG 
equations (O immediately yields the deformed RG equa- 
tions (pTjl . The irrep labels d 2 :(£) are not well-defined for 
£ = 0, enforcing a separate assessment of the contrac- 
tion limit. The reduced BCS Hamiltonian in this limit is 
expressible in terms of hw(l) generators Q17p 

rn m 

H(0) = J2 £ i( rs ° + E (23) 

<=1 i,j=l 

The correspondence between s° and the bosonic number 
operator makes the Hamiltonian essentially a one-body 
bosonic operator 

H(0) = J2 (-' A, + yV^h) + E £ ^' (24) 

i,j=l i—l 

This Hamiltonian can be brought in diagonal form by 
means of a simple unitary transformation of the boson 
operators b\ = Ylk=i ^4 



(25) 



with the cigenmodes fuxii being the roots of the secular 
equation 



E 



9 



2 4 2ei — ftw 

i—l 



0. 



(26) 



This result agrees with the £ — > contraction limit of 
the deformed RG equations (f2"Tj) as lim^o = j^ti, 

pointing out that the su(2)j part (£ G]0, 1]) and the 
hw(l) limit (£ = 0) are consistent within the pseudo 
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deformation formalism. A £ = multiphonon state can 
be represented in multiple ways 

TV m 

WMO)})) = II S l(a)l^> = IK^)"*^) (27) 

a=l fc=l 

fc=l \ i=l 1 V 

with fc(a) the integer-valued function defined as before, 
selecting which eigenmode k with energy hu>k is assigned 
to the ath phonon. The total number of phonons with 
eigenmode huJk in the multiphonon state is given by v/-, so 
we can represent a multiphonon state (|27j) by means of a 
vector (yi, i/2, ... , i/ m ). This definition is consistent with 
the definition given in Sec. IIII Al with the only difference 
that the eigenmodes are here genuine bosons instead of 
generalised pairs. 

In conclusion, the RG equations (J7J of the reduced 
BCS Hamiltonian (TTJ) can be continuously transformed, 
by means of a pseudo deformation of the quasispin al- 
gebra (|13p. into the secular equation ([26]) of a bosonic 
one-body Hamiltonian (|25p bearing a strong resemblance 
with the secular equation (fTU|) of the pp-TDA. The cor- 
respondence between the secular equations (fTO)) and (|26|) 
can be understood from a physical point of view. As the 
eigenmodes of the pp-TDA are one-pair excitations, they 
behave essentially bosonic, so they will be structurally 
equivalent to the eigenmodes of the bosonic Hamiltonian 
([2%]) . The only difference between (JTUJ) and (|26j is the 
presence of the seniority Vi in the pp-TDA. This is a direct 
consequence of how the Pauli exclusion principle is ap- 
plied in both formalisms. Throughout the construction, 
the operator hi was required to retain its interpretation 
as a number operator. Therefore, the number of parti- 
cles occupying the pair vacuum \6) equals the number of 
unpaired particles (Y2iLi v i)i independent from the de- 
formation parameter. The regular fermion pairs will feel 
the presence of unpaired particles, occupying part of the 
free Hilbert space for the fermion pairs, in contrast to the 
pure bosons, who are insensitive to the presence of any 
other particles. 

The advantage of the present construction with the de- 
formed RG equations is that one can make a one-to-one 
connection between the exact eigenstates of the reduced 
BCS Hamiltonian and the multiphonon states by adia- 
batically tuning the Pauli exclusion principle using the 
deformation parameter £. This will be discussed next by 
means of the picket-fence model introduced in the previ- 
ous subsection. 



D. Tuning the Pauli principle 

The deformation parameter £ of the pseudo deformed 
quasispin algebra (|13[) can be used to tune adiabatically 
the effect of the Pauli exclusion principle. So, it is now 
possible to gradually turn off the Pauli principle in a 



given exact eigenstate, and observe unambiguously to 
which multiphonon state it will evolve. The other way 
around, the exact eigenstates can be reconstructed from 
physically relevant multiphonon states by reintroducing 
the exclusion principle. 

Before applying these ideas to the picket-fence model 
of Sec. IIII Al it is instructive to study a few limiting cases. 
Much of our understanding of the general features of the 
RG equations has been derived from these limiting cases, 
such as the weak-coupling [2(| and stron g-co upling 
regime [33(, or the thermodynamic limit [34l. l35j. The 
RG variables in the weak- and strong-coupling regime 
are related to the roots of the associated Laguerre poly- 
nomials using the elegant Hcine-Sticltjes (HS) correspon- 
dence [3(1 [37|- In general, the HS correspondence relates 
polynomial solutions of a type of 2nd-order ordinary dif- 
ferential equations with the roots of a set of non-linear 
coupled algebraic equations. The RG equations can be 
nicely embedded within the HS formalism, which has lead 
to recent developments in the numerics of the RG equa- 
tions (lil. [38l - l40| . For some selected cases, the correspon- 
dence leads to simple approximate expressions for the RG 
variables in terms of the roots of special functions, as will 
be discussed next. Details can be found in Appendix [Bl 

The Heisenberg-Weyl algebraic structure of the con- 
traction limit £ = is quite different from the effective 
quasispin structure for £ g]0,1]. Therefore, it would 
be interesting to see how the RG variables change as 
the Pauli exclusion is turned on infinitcsimally little 
(0 < £ <C 1). Let (v\, v-x-, ■ . ■ , u m ) be the distribution of a 
multiphonon state ([2~T)l in the contraction limit. In this 
limit we start from x a (0) = huJk( a ) for the RG variables, 
with k(a) the integer-valued function as defined before. 
Using the HS correspondence, it is shown in Appendix 
IB II that infinitesimally cranking up the Pauli principle 
(£ <C 1) affects the RG variables as follows 

x a (0 = hw Ha) + iJ^z VKahl{a) + O(0, (29) 

with z v% i the Zth root of the Hermite polynomial H u (z), 
and cifc = i (2e -hio k y 2 a ^-independent constant. The 
integer- valued function 1(a) picks a Hermite root for each 
RG variable, with the only condition that all roots must 
be distinct. In other words, the degeneracy in the Vk 
phonons with eigenmode huJk is immediately lifted by 
turning on the deformation parameter £. The dispersion 
in the RG variables scales with y/^, pointing out that 
the transition from hw(l) into su(2)f, associated with 
the introduction of the Pauli principle, is both analytic 
and non-perturbative. Also, ajt is a positive constant, 
and the roots of the Hermite polynomials are all real- 
valued, showing that all (but one) RG variables immedi- 
ately become complex valued for even (odd). There- 
fore, we can conclude that it is in fact the Pauli princi- 
ple pushing those pairs associated with the same eigen- 
mode into the complex plane. It is worth mentioning 
that a HS correspondence with the Hermite polynomi- 
als has been identified earlier in the weak-coupling limit 
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of the ID 5-interacting Licb-Liniger bose gas [41| . More- 
over, a connection with the strong-coupling regime of the 
reduced BCS Hamiltonian has been established how- 
ever there appears to be no connection with the pseudo 
deformation scheme of the present article. 

The previous analysis demonstrates the behavior of the 
RG variables as they move away from the bosonic con- 
traction limit. However, the description (|29[) is only valid 
in the vicinity of the contraction limit, and it is impos- 
sible in general to extend this to the fully restored qua- 
sispin limit £ = 1. Nonetheless, it is known that the 
RG variables can be related to the roots of the associ- 
ated Laguerre polynomials in both the weak- and strong- 
coupling limit [111 HH. It will be demonstrated in Ap- 
pendix IB 21 that this result is extendible along the full 
path £ € [0, 1], allowing for a 1-to-l adiabatic connection 
between the exact ground state and its bosonic counter- 
part in the contraction limit. The RG variables are given 
by 



x a (0 = -g£yi( a ) + 



(30) 



with yi shorthand for the Ith root of the associated La- 
guerre polynomial i^ 2d(C)_1 (y), and d(£) = YT=i d i(0- 
The value — 2<i(£) — 1 is negative, so all roots are non-real, 
except for at most one real root when N is odd [13, S3 . 
The description for the RG variables (|3T)|) is in principle 
valid for any finite value of £. Caution is required for the 
contraction limit £ — > 0, as the parameter — 2c£(£) — 1 of 
the associated Laguerre polynomial becomes singular. It 
can be shown that the roots of the associated Laguerre 
polynomial L^(y) with |fc| — > oo behave approximately 
like H El 



Vl 



2kz n i, 



(31) 



with z n .i the Zth root of the Hermite polynomial H n (z). 
Using this relation in (l30l . combined with the fact that 



IQ. 



£,hvi is well-defined and finite, we ob- 



tain the near-contraction limit for the ground state of 
the strong-coupling regime 



M)=9 



£flz N 



1(a) 



(32) 

with f2 = YjiLi This result is consistent with (|2"5|) . as 
the connecting multiphonon distribution of the ground 
state is given by (N, 0, . . . , 0), with lowest (collective) 
eigenmodc 



a 



(33) 



The analysis of the strong-coupling regime within the HS 
formalism shows that the exact ground state of the re- 
duced BCS Hamiltonian can be unambiguously related 
to the (TV, 0, . . . , 0) collective multiphonon state by adia- 
batically switching off the Pauli exclusion principle via a 
pseudo deformation of the quasi spin. It should be men- 
tioned that the HS correspondence allows for a qualita- 
tive description, which is only correct within the validity 



domain of the approximation. For an exact assessment 
of how the exact exact states are related to the bosonic 
counterpart, it is necessary to numerically solve the de- 
formed RG equations (|2"Tj) . This was done for the picket- 
fence model [20|, described in the previous section, for 
a selected set of coupling constants g. The results are 
presented in Fig. 0] and Fig. [5] Fig. 0] shows the path 
of the RG variables in the complex plane as the pseudo 
deformation £ is tuned down from 1 to 0. This is done 
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FIG. 4. The path of the deformed RG variables {x(£)} in the 
complex plain for the picket-fence model of Sec. IIII Al with 
\g\ = 0.650 (a), and |g| = 5.000 (b). the path, given by full 
lines for each variable, starts from the exact RG variables 
{a;(l)} (filled circles), and ends at to the bosonic eigenmodes 
{:e(0)} = {huj} (open circles). The near-contraction approxi- 
mation is given by dashed lines, and the strong-coupling limit 
approximation is given in dotted lines (only for (b)). The 
inset in both figure shows a close-up of the path around 
the near-contraction limit. It can be remarked that the 
near-contraction limit predictions and the numerically ob- 
tained deformed RG variables nearly coincide at the near- 
contraction limit for \g\ = 5.000. The single-particle levels 
at 2ei (i = l..m) are also included in dashed-dotted lines as 
reference points for the multiphonon states. 

for a system in the intermediate- (\g\ = 0.650 in Fig.@Ji) 
and strong-coupling (\g\ = 5.000 in Fig. 0Jd) regime. It 
is readily observed that the RG variables evolve towards 
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real and degenerate points in the complex plane, associ- 
ated with the multiphonon states. For the intermediate 
case, the associated multiphonon state is (1, 1, 4, 0, ... ), 
whereas the strong-coupling case connects with the col- 
lective (6, 0,0,...) multiphonon state. In both cases, the 
exact path has been confronted with the predictions of 
the near-contraction limit, up to 0(£) (see Appendix IB 11 
on how the approximate values of Eq. (|29[) can be im- 
proved up to C(£)). It can be concluded from the insets 
in Fig. H] that the near-contraction limit of the deformed 
RG equations is well understood using the HS correspon- 
dence. In addition, the results for the strong-coupling 
case in Fig. [4j} are compared with the approximate ex- 
pression (|30|) in the strong-coupling limit. The general 
trend of the RG equations in the complex plane is rather 
well reproduced, taking into account that the errors are 
expected to be of order 0(g~ 2 ) with respect to the or- 
der O(g) of the RG variables. So, it is confirmed that 
the RG variables in the strong-coupling regime indeed 
evolve towards a fully collective multiphonon state, as 
anticipated. 

It can be inferred from Fig. |4] that the exact ground 
state of the reduced BCS Hamiltonian can connect to dif- 
ferent multiphonon states for different coupling constants 
g. This has been investigated in detail for the ground 
state of the picket fence model of section [Til Al where the 
interaction strength has been varied from g = —2.000 to 
g = 0.000. The results of this study can be found in 
Figures [5] & [TJ It was observed that the exact ground 
state in the strong-coupling regime is always connected 
with the (6,0,..., 0) multiphonon state (see Fig. 2] and 
Fig. [3 for respectively \g\ = 5.000 and \g\ = 1.000), 
until the critical interaction strength \g c \ = 0.890 was 
reached as a lower limit. Crossing this critical point gave 
rise to a change in multiphonon character in the con- 
traction limit. Instead of the fully collective state, the 
exact state chose to connect with the (5, 1, 0, . . . , 0) mul- 
tiphonon state (see Fig. ET) for \g\ = 0.850). This can 
be understood by realizing that one of the complex con- 
jugate pair of RG variables changes into two real-valued 
variables at the critical interaction strength. In order to 
show the subtleties involved with the change of the mul- 
tiphonon state around the critical interaction strength 
g, the path towards the contraction limit was plotted in 
Fig. Eg and Fig. Eh for \g\ = 0.890 and \g\ = 0.891 re- 
spectively. The resolution of Fig. El is hardly sufficient to 
resolve the two critical RG variables in both the figures, 
however the path to the multiphonon states is clearly dif- 
ferent. Remarkably, the next g where the multiphonon 
state changes character is not at the next critical interac- 
tion strength \g c \ = 0.603. It was observed that the next 
two instances occurred at \g\ = 0.845 and \g\ = 0.719 
where the multiphone state changes to (1, 5, 0, . . . , 0) and 
(1, 1, 4, 0, . . . , 0) respectively, which can be interpreted 
as a simple rearrangement of the two real-valued and 
four complex-valued RG variables into a less collective 
configuration. The path towards the (1, 5, 0, . . . , 0) and 
(1, 1, 4, 0, . . . , 0) multiphonon states are illustrated for 



|S| = 0.800 and \g\ = 0.650 in Fig. EH and Fig. Et- 
The near-contraction limit of the latter figure has al- 
ready been discussed in detail in Fig. |4ji. The next 
point where the multiphonon configuration alters is at 
the critical point \g c \ = 0.603, where the (1, 1, 4, 0, . . . , 0) 
changes into (1, 1, 1, 3, 0, ... , 0), which can again be re- 
lated to the breaking of a complex conjugate pair of 
RG variables into two real- valued RG variables. The 
path to the (1, 1, 1, 3, 0, . . . , 0) multiphonon state is il- 
lustrated for \g\ = 0.600 in Fig. Et- Next, there is again 
a change in multiphonon distribution related to a rear- 
rangement of the real- and complex- valued RG variables 
at \g\ = 0.580 to a less collective (1, 1, 1, 1, 2, . . . , 0) (see 
Fig. Eb)- Finally, the last change occurs at the critical 
value \g c \ = 0.393, below which each exact state cor- 
responds to a (1, 1, 1, 1,1,1,0,..., 0) multiphonon state 
(as illustrated for \g\ = 0.200 in Fig. EM- The phys- 
ical interpretation of this configuration is a simple fill- 
ing of the single-particle levels with pairs until the Fermi 
level is filled. The situation of the picket-fence model 
of subsection IIII Al is summarized in Fig. [JJ and Table Q] 
The coupling constants for which the connecting multi- 
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TABLE I. The coupling strengths \g\ where the connect- 
ing multiphonon character changes for the picket-fence 
model of Sec. IIII Al The given multiphonon distributions 
(yi, i/2, ■ ■ ■ , Vm) in the row of each \g\ is the distribution into 
which the change occurs for an increasing \g\. The label r 
and c denote whether the coupling constant corresponds to 
respectively a rearrangement or a critical RG variable. These 
interaction strengths are visualized in Fig. [TJ 

phonon distribution changes have been denoted by dot- 
ted and dashed lines in Fig. [TJ with the distinction that 
the dashed lines denote the critical interaction strengths 
associated with a recombination of RG variables into 
complex conjugate pairs, whereas the dotted lines indi- 
cate the values of g where the connecting multiphonon 
changes without such a recombination. The exact values 
where the multiphonon distributions change are given in 
Table U up to three decimal places. 

Now, it would be interesting to put the present clas- 
sification of the ground state into different multiphonon 
distributions against the overlaps of the exact ground 
state with the pp-TDA multiphonon states (ip{x}\ip{fiLu}) 
of Sec. IIII Al (see Fig. The present discussion an- 
swers the question of Sec. IIII Al whether it is possible 
to assign a change of multiphonon character to the crit- 
ical point \g c \ = 0.603. However, caution is required 
identifying the contracted multiphonon states (|27|) with 
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FIG. 5. The path of the deformed RG variables {x(£)} in the complex plain for the picket-fence model of Sec. MI Al with 
\g\ = 0.200 (a), \g\ = 0.500 (b), \g\ = 0.600 (c), \g\ = 0.650 (d), \g\ = 0.800 (e), \g\ = 0.850 (f), \g\ = 0.890 (g), \g\ = 0.891 (h), 
and \g\ = 1.000 (i). The path is given starting from the exact RG variables {ie(1)} (filled circles) and ending at the bosonic 
eigenmodes {:r(0)} = {fru}} (open circles). The connecting multiphonon distribution (yi, V2, ■ ■ ■ , v m ) is included in each figure, 
and the single-particle levels at 2si (i = 1 . . . m) are also given in dashed-dotted lines as reference points for the multiphonon 
states. 



the pp-TDA multiphonon states (fTTj) . For the sake of 
comparison, the overlaps of the TDA phonon states 
with the exact ground state for the distributions given 
in Table U have been plotted as well in Fig. [2] It 
can be seen that the identification (contraction-limit vs. 
pp-TDA multiphonon states) works rather well in the 
weak- and strong-coupling regime. Indeed, for what the 
weak-coupling regime is concerned, the domains where 
the exact ground state connects with the multiphonon 
states (1,1,1,1,1,1,0, ...,0) and (1, 1, 1, 1, 2, 0, . . . , 0) 
correspond largely with the domains where the 
(1, 1, 1, 1,1,1,0,..., 0) and (1, 1, 1,1,2,0,..., 0) pp-TDA 
multiphonon states give the maximum overlap (compare 
Table U with Fig. [2]). Similarly, in the strong-coupling 
regime, there is a correspondence between the domains 
of the connecting multiphonon states (5, 1, 0, . . . , 0) and 
(6, 0, ... , 0) with respectively the domains where the pp- 
TDA states with eigcnmode distribution (5, 1, 0, . . . , 0) 
and (6,0, ...,0) give the maximum overlap. However, 
the connection is not unambiguous in the intermediate 



regime. Closer inspection of Fig. [5] reveals that the over- 
laps of the pp-TDA states with eigcnmode distribution 
(1, 1,1, 3,0,. ..,0), (1, 1,4,0,. ..,0) and (1, 5, 0, . . . , 0) 
(c/r. Table IJ) may grow relatively large in the intermedi- 
ate regime, but they never result in the maximum over- 
lap for a given interaction strength. This points again 
towards the intricate role the Pauli principle plays in the 
intermediate regime. The exclusion principle is fully ac- 
tive in the calculation of the overlaps of the exact with 
the pp-TDA multiphonon states, because both states are 
living in the same Hilbert space. So, the TDA state with 
maximum overlap will be the result of a subtle trade off 
of constructing a highly collective multiphonon state (i.e. 
more phonons with eigcnmode hull) against the action of 
the Pauli principle, washing away part of the collectivity. 
Within the framework of the pseudo-deformed quasispin 
algebra, this trade-off is irrelevant as the Pauli principle 
is totally suppressed in the contraction limit. 

The critical values g c are clearly related to a change 
of multiphonon character in the contraction limit. How- 
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ever, as is clear from the previous discussion, in addition 
to these critical values, a few other values of g have been 
found where the connecting multiphonon state changes 
character. It would be interesting to investigate a possi- 
ble relation between these changes in multiphonon char- 
acter and higher-order singular points in the RG vari- 
ables. For this purpose, the derivatives of the RG vari- 
ables have been plotted in Fig. [51 with the antici- 
pation that some additional structure could be found at 
these particular values of g. Unfortunately, as can be 
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exact RG states ((6|) with the corresponding contracted 
multiphonon state ([27]). A valuable question would be 
how the connection works in the opposite direction, i.e. 
given a multiphonon state (|27[) . is there a physical state 
([6]) , corresponding to an exact eigenstate of the reduced 
BCS Hamiltonian ([lj? The Hilbert space H^ = o spanned 
by all possible multiphonon states is much larger than 
the Hilbert space T-L^—i spanned by all possible hard-core 
bosonic states. For instance, the dimensions of a system 
of N pairs living in m twofold degenerate (fl; = 2) single- 
particle levels are given by 

dimH^o = ( m + ft 1 J , dimH ?= i = ( ^ V 

(34) 

Therefore, we may expect that there are multiphonon 
states for which there is no corresponding exact state. 
Indeed, as an example, adiabatically turning on the Pauli 
principle starting from the collective (6,0, ... ,0) multi- 
phonon state did not result into an exact eigenstate at 
£ = 1 for \g\ < 0.891. The solution broke down at a 
given £ < 1 along the path. This is somewhat unfortu- 
nate, because the quasi-deformed algebra approach could 
have the potential to act as a convenient RG solver, com- 
plementary to other approaches @, Q HI, [H S3] . Mul- 
tiphonon states are considerably more tractable from a 
numerical point of view than the RG equations, and the 
Pauli principle can be reintroduced in a very controlled 
fashion. The fact that the Hilbert space of the genuine 
bosons is larger than the Hilbert space of the fermionic 
pairs, would imply that many multiphonon states may be 
tried before the corresponding exact state is found. Nev- 
ertheless, as has been shown in the present manuscript, 
as well as in previous studies 0, [47|, physical insight 
and intuition can prove very helpful in finding the mul- 
tiphonon states corresponding to the exact ground state 
and excited states [2(| Esl • As a matter of fact, the ex- 
act ground states in the present work have been obtained 
using this idea of reconstruction the exact ground state 
from the connection multiphonon state. 



FIG. 6. The real part (a), imaginary part (b), and the total 
sum (c) of the derivatives of the exact RG variables with re- 
spect to the interaction strength g. The singular interaction 
strengths are highlighted by dotted and dashed lines, where 
dashed lines denote critical values \g c \ and dotted lines high- 
light the values of \g\ where a multiphonon rearrangement 
occurs (see Fig. [TJ and Table [I}. 

seen from the figure, no hints of a possible higher-order 
singular point could be found at the mentioned values of 
g. Therefore, the nature of these additional points re- 
mains largely unclear. It is worth mentioning that none 
of these points were found in previous a naly sis' within the 
pseudo-deformed quasispin formalism (46ll47|. including 
a study of the numerical test-case model of Rombouts et. 
al. |3 • 

The discussion has focused thusfar on connecting the 



IV. CONCLUSIONS 

In conclusion, an adiabatic connection has been made 
between the hard-core bosonic eigenstates of the reduced 
BCS Hamiltonian and the multiphonon states charac- 
teristic to bosonic approximations, using a pseudo de- 
formation of the su(2) quasispin algebra. The pseudo 
deformation is not a genuine deformation because it 
can be transformed into an effective quasi-spin algebra 
with rescaled irrep labels. Only in the contraction limit 
of the deformation parameter transforms the deformed 
su(2) into a Heisenberg-Weyl algebra hw(l), the fun- 
damental spectrum generating algebra of bosonic exci- 
tations. It is shown that the reduced BCS Hamilto- 
nian remains Richardson- Gaudin integrablc along the de- 
formation path of the pseudo-deformed quasispin alge- 
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bra with a corresponding set of deformed Richardson- 
Gaudin equations. This set of deformed Richardson- 
Gaudin equations become equivalent to the seniority-free 
secular equation of the pp-Tamm-Dancoff Approxima- 
tion. The near-contraction and strong-coupling limit of 
the deformed Richardson-Gaudin equations have been in- 
vestigated using a Heinc-Stieltjcs correspondence and an 
example using a so-called picket-fence model is discussed 
in detail. The present approach allowed for an interpreta- 
tion of the critical interaction strengths g c as those points 
where the connecting multiphonon states change collec- 
tive nature. Also, the potential of the adiabatic con- 
nection as a Richardson-Gaudin solver has been briefly 
discussed, complementary to recent developments in the 
numerics of the Richardson-Gaudin equations. Finally, 
it would be interesting to investigate how this pseudo- 
deformation scheme ispositioned within the huge realm 
of bosonic mappings [48( and bosonic approximations to 
pairing dynamics of fermionic systems [49| . 
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Appendix A: Deriving the RG equations 

It is straightforward to derive the RG equations using a 
commutator scheme @, [26| . In this scheme, the Hamilto- 
nian (TTJ) (or Eq. (fT9|) in the deformed case) is commuted 
through the product state ((6|) (or Eq. ([20]) ) until it acts 
on the pair vacuum. The RG equations are obtained 
by requiring the non-diagonal by-products to vanish. It 
is convenient to introduce the following notation. Let 
Oi be a level- dependent operator with the Roman index 
i = 1 . . . m denoting the level. The associated Gaudin 
operator is given by 



On 



1=1 



Oi 



(Al) 



with the Greek label a shorthand for the RG variable x a . 
Further, it is also convenient for notational purposes to 
introduce the operator O = X)"=i Oi- Using this nota- 
tion, the Schrodinger equation becomes 



JV 



JV 



h{o n sum = m n Qim- ( A2 ) 



a=l 



a=l 



The general deformed case (as defined in Sec. IIII B| will 
be used in the present Appendix, as the regular RG for- 
malism can be obtained from the £ = 1 limit. By com- 
muting the Hamiltonian through the product state, we 
derive the following expression 

N N 

H(0 [J St(0\6) = I] Sl(0H(0\e) (A3) 

a— 1 a— 1 

N N 

JV JV N 

+E E II si(oimo,sl(0},s\(0¥), 

0=1 7=0+1 a#/3,7 

where the property that the double commutator with the 
Hamiltonian (|19l) commutes with all creation operators 
has been used in the last line. The single and double 
commutators can now be calculated explicitly using the 
deformed algebra p^|) . 



+ SH0 [i-ff(26$2(0 + (£-i)3 n /s 



5^(0- 5t (0 



(A4) 



(A5) 



The action of these commutators on the pair vacuum 
state simplifies Eq. (|A3I) to 



m) n & a (m 



a=l 



N 



0=1 



N 

E 

0=1 



i=l 



I sum 



l + 2# ? (0-2 S ^ 



a=l 



1 



7#/3 



N 



(A6) 



taking into account that 2£d a + (£ — l)gfi Q = 2£g? q (£) 
(see Eq. (fT()]) ). The crucial observation is now that the 
product state [] , S^\6) will only be an eigenstate iff 
the non-diagonal terms vanish identically. This is the 
case if 



N 



1 + 2 ff ^(o -2 5 eE 



1 



= 0, 



7/0 



V/3 = 1 . . . N, 
(A7) 
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which constitutes identically the set of deformed RG 
equations (|21[) . It is worth pointing out that the present 
general procedure logically applies to the regular RG 
formalism for the BCS Hamiltonian (£ = 1) as to the 
bosonic contraction limit (£ = 0). Moreover, it has been 
used to derive the Bethe Ansatz of the Dicke model, 
describing a mixture of hard-core bosonic and genuinely 
bosonic pairs. 



Appendix B: Limits of the deformed RG equations 

In some limiting cases, the solutions of the RG equa- 
tions can be related to the roots of special functions, us- 
ing the Heine-Stieltjes correspondence [H, H3 • This will 
be discussed in more detail in the present Appendix for 
the near-contraction limit and the strong-coupling regime 
over the whole regime of the deformation parameter. 



1. The near-contraction limit 

Assume we have a multiphonon state (|27[) in the 
contraction limit, with a given eigenmode distribution 
(yi, i>2, ■ ■ ■ , I'm)- The RG variables in this limit arc given 
by x a (0) = huj k r a \, with hu k the solutions of the sec- 
ular equation (|26p and k(a) an integer-valued function 
picking the fcth eigenmode for the ath pair. Making use 
of (fT5|) . the deformed RG equations (|2"Tj) can be written 
explicitly as 



Now, £ can be chosen small enough for a meaningful series 
expansion up to first order in y 7 ^. 



1- 



fe(a) 



(B4) 
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1 m 



2 frj- (2 £i 

2—1 V 

O(0 w0, (Va=l...iV), 



E 



.a) 



X /3 x 



(1) 



with the summation over (3 limited to those pairs for 
which k((3) = k(a). The key observation is that the terms 
in 0(1) vanish by definition because twj k has been chosen 
as such that it resolves the bosonic secular equatio n (1261) . 
As a result, we obtain the Stieltjes equations [43l l5ll| 



Ofc(a) 



,(1) - 



E 



(i) 



Va = 1 . . . N 



(B5) 



with ai = 5 X)i=i (2e -huii) i 1 Because the summation /3 
runs only over the RG variables with k{0) — k(a), the full 
TV-pair problem breaks down into separate blocks of RG 
variables with equal x a = fvjj k i a \, for which the Stieltjes 
equations can be solved independently from each other. 
This can be done using the Heine-Stieltjes connection 
[43l l5ll ] , of which the main ideas are summarized below. 
Assume we want to solve the Stieltjes equations for v 
variables x^ 1 (a = Let's introduce the Heine- 

Stieltjes polynomial, 

P(x) 



(B6) 



N 



m c In 



. 2 £i - Xot Q . XR Xa. 

8=1 P^a r 



Assume that the RG variables in the near-contraction 
limit are given by 

x a (0 = hu k(a) + + O(0- (B2) 

We need to make a distinction in the summation over /3 
in equation (|B1|) between the pairs xp{£) with the same 
eigenmode as x a (£) (k((3) = k(a)) and the pairs 2/3'(£) 
with a different eigenmode (k((3') ^ k(a)) 



wE 



—[ 2e l - hhJ k{a) - y/lXa 1 ' 

-2tf E 



(B3) 



2 9^l E (i) ~ (1) (Va=l.JV). 



k(p)=k(u) 



defined as the polynomial of degree v with the solution 
of (|B5P as the roots. It is straightforward to derive that 



P"(xP) v 1_ 



(Va = 1 . . . v) (B7) 



As a result, the Stieltjes equations (|B5|) can be rewritten 

as 

«fc4 1) ^'(4 1) )+ p "(4 1) ) = I (Va=l...i/) (B8) 

Now, we can define the polynomial 

Q{x) = a k xP'{x)+P"{x) 1 (B9) 

which has the same degree v as P(x), and the same set 
of roots (see equation (|B8jl ). so Q(x) must be a multiple 
of P{x). It is straightforward to see from the definition 
of Q{x) that the multiplication factor is a k v. Therefore, 
we obtain the equation 



a k P'(x) + P"{x) = a k vP(x). 



(BIO) 



A change of variable z = —i^Ja k /2x transforms this 
equation into the Hermite differential equation [52I ] 



H"(z) - 2zH'(z) + 2vH(z) = 0. 



(Bll) 
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In conclusion, the solution of the Stieltjcs equations (|B5[) 
can be related to the roots {z„} of the Hermite poly- 
nomials H v (z). Consequently, the RG variables in the 
near-contraction limit are given by 



x a (0 = hu k (a) +i 



2g 



0(0, (B12) 



with z v> i the Ith root of the Hermite polynomial H„(z). 

It is understood from equation (|B12I) that the first- 
order correction in is purely imaginary. Therefore, if 
a qualitative description of the real part of the RG vari- 
ables near the contraction limit beyond 0(1) is desired, 
it is necessary to include higher-order terms in the ap- 
proximate solution of the deformed RG equations ([Bljl 



x a (0 = hui k{a) + y/ix& + (B13) 

In addition to the Stieltjcs equations at order 0(\/£), we 
obtain a linear set of equations for the correction vari- 
ables {x (2) } at order 0(£) 



afc(a)a4 2) 



f3^a ( x $ x a ) 



E 



KP)J=k(a) 
-&fc(a)(z (1) ) 2 +c fc(Q ) +B(k(a)), 



(B14) 



with ai defined previously, b t = \ YaL\ (2 Ei -L,) j ' and 

with the summation ft' restricted to the pairs for which 
k(ft') 7^ k(a), is a functional of the integer- valued func- 
tion k(a). It can be verified, taking the Stieltjcs relation 
(|B5jl for the variables {x^ 1 )} into account, that the solu- 
tion to this set is given by the simple relation 



,(2) 



— y x *) +^72 — K(a)-l) 
6a k{a) 6a k ( a ) 

c k(a ) B(k{a)) 



(B15) 



In conclusion, the 0(£) corrections are purely real, and 
are directly related to the 0( v / £) corrections. These cor- 
rections have been taken into account for the comparison 
with exact results in Fig. 



2. The strong-coupling limit 

A similar analysis using the HS correspondence can be 
performed for the strong-coupling regime (2(1 [H| . The 
strong-coupling case of the present subsection is different 
with respect to the near-contraction case of the previous 
subsection as the present analysis is valid for any value 
of the pseudo deformation parameter £ , whereas the pre- 
vious analysis was only legitimate near the contraction 



limit. It is instructive to start with the 1-pair problem, 
with the deformed RG equation (|2"Tj) given by 



0. 



(B16) 



The ground-state solution of this equation in the \g\ — > oo 
limit is 



x a (0 = 2^d(0 



(B17) 



with = YliLi This result suggests that the 

ground state of the iV-pair problem will be described by 
RG variables of the following form (the £ dependency will 
be omitted for notational simplicity) 

x a =rf+x«+0(±), Va = l...JV. (B18) 

Inserting this in the deformed RG equations (1211) . we 
obtain 



1 + 29J2 



— 1 fj^a %a 

N 



E 



i 



X 1 ) 



(B19) 
0. (B20) 



A Taylor expansion up to next-to-leading order simplifies 
this equation to 



' jo) 



N 



* E 



i 



p=l,jta x fi 



(0) _ _(0) 



(B21) 



(i) 



(xL 0) ) 2 



r (i) 

x^ 



(xL 0) ) 2 



£ E 



Xct ) 



+ 0(4 r )=O,Va. 



We require each order to vanish. First, it can be verified 
that 



Va = 1 . . . N 



(B22) 



is a solution for the 0(^) term in equation (|B21[) . in- 
dependent from the solution {x' - 1 } of the 0(1) term. 
Second, the solution {x^ '} of the 0(1) term can be 
found using the HS correspondence, similar to the near- 
contraction limit (Appendix IB 1[) . Again, by introducing 
the HS polynomial with the values {x^ ^} as roots 



N 



(B23) 



it is possible to rewrite the equation resolving the 0(1) 
term as 

x^P'(x^)-2^d(0P'(x^) + ^ ) P''(x^) = 0. (B24) 
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The polynomial 

Q(x) = xP'(x) - 2£d(£)P'(x) + £xP"(x), (B25) 

has the same degree N and roots {x^} as P(x). There- 
fore it must be a multiple of P(x), with the multiplication 
factor easily determined as N 

xP'(x) - 2£d(£)P'(x) + &P"(x) = NP(x). (B26) 

A change of variables x = — £?/ transforms this equation 
into the differential equation defining the associated La- 
guerre polynomials L^{y) of degree k = — 2d(£) — 1 [52[ 



y f3^1 + (fc + 1 - v) d JMyl + NLk {y) = o. (B27) 

dy dy 

As a result, the variables {x^} can be related to the 
roots of the associated Laguerre polynomials L^(y) of 
degree k = —2d(£) — 1, and the solution of the deformed 
RG equations in the strong-coupling limit can be written 
as 

sa(fl - -gftftfr) + Er =^g €)2Ei +0(|), (B28) 

with j/; shorthand for the Zth root of the associated La- 
guerre polynomial L^(y) of degree k — — 2d(£) — 1. 
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